Traffic of single-headed motor proteins KIFIA: eff'ects of lane changing 
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KlFlA kinesins are single-headed motor proteins which move on cylindrical nano-tubes called 
microtubules (MT). A normal MT consists of 13 protofilaments on which the equispaced motor 
binding sites form a periodic array. The collective movement of the kinesins on a MT is, therefore, 
analogous to vehicular traffic on multi-lane highways where each protofLlament is the analogue of 
a single lane. Does lane-changing increase or decrease the motor flux per lane? We address this 
fundamental question here by appropriately extending a recent model [Phys. Rev. E 75, 041905 
(2007)]. By carrying out analytical calculations and computer simulations of this extended model, 
we predict that the flux per lane can increase or decrease with the increasing rate of lane changing, 
depending on the concentrations of motors and the rate of hydrolysis of ATP, the "fuel" molecules. 
Our predictions can be tested, in principle, by carrying out in-vitro experiments with fluorescently 
labelled KIFIA molecules. 



\ Members of the kinesin superfamily of motor proteins 
move along microtubules (MTs) which are cylindrical 
nano-tubes [H, [2]. A normal MT consists of 13 protofil- 
aments each of which is formed by the head-to-tail se- 
quential lining up of basic subunits. Each subunit of a 
protofilament is a 8 nm heterodimer of a-j3 tubulins and 
provides a specific binding site for a single head of a ki- 

] nesin motor. Often many kinesins move simultaneously 
along a given MT; because of close similarities with vehic- 
ular traffic 3 , the collective movement of the molecular 
motors on a MT is sometimes referred to as molecular 
motor traffic [1, i, 0, 0, l . 

The effects of lane changing on the flow properties of 
vehicular traffic has been investigated extensively using 
particle-hopping models [s*] which are, essentially, appro- 
' priate extensions of the totally asymmetric simple exclu- 
sion process (TASEP) 0, O, |il|. Models of multi-lane 
TASEP, where the particles can occasionally change lane, 
have also been investigated analytically [l^, E^l- Two- 
lane generalizations of generic models of cytoskeletal 
molecular motor traffic have also been reported [H, [lB| • 

Recently a, quantitative theoretical model has been de- 
veloped [la [i3| (from now onwards, we shall refer to it 
' as the NOSC model) for the traffic of KIFIA proteins, 
which are single-headed kinesins [lM[S[13|i by explicitly 
\ capturing the essential features of the mechano-chemical 
I cycle of each individual KIFIA motor, in addition to their 
steric interactions. In this communication we extend the 
NOSC model by adding to the master equation all those 
terms which correspond to lane changing. Solving these 
equations analytically, we address a fundamental ques- 
tion: does lane changing increase or decrease flux per 
lane? We show that the answer to this question depends 
on the parameter regime of our model. We establish the 
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levels of accuracy of our analytical results by compar- 
ing with the corresponding numerical data obtained from 
computer simulations of the model. We interpret the re- 
sults physically and suggest experiments for testing our 
theoretical predictions. 




FIG. 1: (Color online) A schematic representation of our 
model, together will the allowed transitions and the corre- 
sponding rate constants. 

The equispaced binding sites for KlFlA on a given 
protofilament of the MT are labelled by the integer in- 
dex z (z = 1, ...,L). We use the integer index j to label 
the protofilaments; the position of each binding site is 
denoted completely by the pair Because of the 

tubular geometry of the MT, periodic boundary condi- 
tions along the j'-direction would be a natural choice. 
We impose periodic boundary conditions also along the 
i-direction, as it not only simplifies our analytical calcu- 
lations but is also adequate to answer the fundamental 
questions which we address in this communication. 

In each mechano-chemical cycle a KIFIA motor hy- 
drolyzes one molecule of adenosine triphosphate (ATP) 
which supplies the mechanical energy required for its 



2 



movement. The experimental results on KIFIA motors 
[Tsl . [TqI . [20j indicate that a simplified description of its 
mechano-chemical cycle in terms of a 2-state model 
would be sufficient to understand their traffic on a MT. 
In the two "chemical" states labelled by the symbols S 
and W the motor is, respectively, strongly and weakly 
bound to the MT. 

In the NOSC model, a KIFlA molecule is allowed to 
attach to (and detach from) a site with rates uja (and 
uj). The rate constant lO(, corresponds to the unbiased 
Brownian motion of the motor in the state W. The rate 
constant cuh is associated with the process driven by ATP 
hydrolysis which causes the transition of the motor from 
the state S to the state W. The rate constants ujf and 
Us, together, capture the Brownian ratchet mechanism 
|2ll . |23 | of a KIFlA motor. Moreover, any movement 
of the motor under these rules is, finally, implemented 
only if the target site is not already occupied by another 
motor. 



The rules of time evolution in the extended NOSC 
model proposed here are identical to those in the NOSC 
model, except for the following additional lane-changing 
rules (see figH]): 

a motor weakly-bound (i.e., in state W) to the binding 
site i on the protofilament j is allowed to move to the 
positions + 1) and (i, j — 1) 

(i) without simultaneous change in its chemical state, 
both the corresponding rates being uju ; 

(ii) with simultaneous transition to the chemical state S, 
the corresponding rate constants being iVfi^ and 
respectively. 

Let Si{j,t) and Wi{j,t) denote the probabilities for a 
motor to be in the "chemical" states S and W, respec- 
tively, at site i on the protofilament j. In the extended 
NOSC model, under mean-field approximation, the mas- 
ter equations for the probabilities Si{j, t) and Wi{j, t) are 
given by 



dS.,[j,t) 
dt 



ua[l - S,(3, t) - W^{j, t)] - L0hS,{3, t) - LUdS,{j, t) -f u;M^{i: t) + C^/W^^-l(j, t)[l - 5,(j, t) - W,{j, t)] 



+LUfi+mj ~ l,t)][l - 5,(j) - Wdj)]+LUfi4W,{j + l,t)][l - S,ij) - W,{j)], 



(1) 



dW,ij,t) 
dt 



i{j,t)-W,+i{j,t)] 
W,-i{j,t)] 



UhS^U, t) - LOsW^ij, t) - LUfW^ij, t)[l - S, 

-iotW.U, t)[2 ~ S,+i{j, t) - W^+i{j, t) - 5,_i(i, t) 
+a;6[I¥.-i(j,<) + I^,+i(j,t)][l - S,{j,t) ~ W,{j,t)] 
+uum{3 - 1, i) + W,(j + 1, t)][l - S,{j,t) - W,{j, t)] 
-u;Mj,t)[2 - 5,(j + l,t)- W,{j + - S,{j ~ l,t) 



W,{3-l,t)] 
ujfi.WdjM^-S,{j~l,t) 



W,{j-l,t)]. (2) 



Rate constant 


numerical value/range (s ^) 




0.1 - 10.0 


<^d 


0.1 




- 250 


LOs 


145 


Uf 


55 



TABLE I: Rate constants of the NOSC model which have 
been extracted from empirical data on single KIFlA experi- 
ments IitIi. 



~ _ rife + + {n, + i)K - \/D 

2K ' ^ ' 

where K — LOd/uja, = '-^h/^f, — ^sl^f-, with 

D ^ m^Kil + K) + {Qh + ns + [Qs - l)K)^- (5) 

The average total density of the motors attached to each 
filament of the MT in the steady state is given by 



In the steady state under periodic boundary conditions, 
S — Si{j,t) and W — Wi{j,t), independent of t and 
irrespective of i and j; from eqs.((T|) and we get 
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2K{1 + K) 



(3) 



~ ^ ~ n,, + + {fis + i)k-Vd + 2 
p = s + w = . (6) 

Using the expressions ([3]) and ([4]) for 5* and W , respec- 
tively, in the expression 

J ^ujfW{l- S ~W) (7) 

for the fiux of KIFlA motors per lane of the MT highway. 
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FIG. 2: (Color online) Flux per lane (and, in the inset, av- 
erage density of the motors on each lane) are plotted against 
ujfi/ujf for a few values of ujs/ujf. Our mean-field predictions 
(labelled MF) are plotted by lines while the discrete data 
points (labelled Sim) have been obtained from our computer 
simulations of the model. 
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FIG. 3: (Color online) Same as in Fig[2] except that the data 
are plotted against K for a few values oi ujfi/ujf. 



we get 



^h + {l + K)^s 



AKil + K) 



(8) 



For graphical presentation of our main results, we use 
the estimates of the rate constants, listed in tablelH which 
were extracted earlier [l^l from empirical data on single 
KIFIA (we use loh = 125 s~^). Since no estimate of 



and LOfi_ are available, we use ^fi+ 



Jfi- 



uifi and 



vary the single parameter ujfi/ujf over a wide range to 
explore the consequences of different rates of lane chang- 
ing. The flux per lane (obtained from ([5])) and the aver- 
age density p (given by ^) are plotted against tofi/ujf in 
Fig. [21 for several different values oiojs/uif and compared 
with the corresponding simulation data. 



FIG. 4: (Color online) Same as in Fig[2l except that the data 
are plotted against uJh/u)f for a few values of K. 



Recall that flux is essentially an average of the product 
of the density and speed of the motors. For sufflciently 
high LUs/ujf, the density p is small even in the absence 
of lane changing {lu fi — 0) and, consequently, the motors 
feel hardly any steric hindrance; increasing tofi/uf in this 
regime of tUs/uif has very little effect on the average speed 
of the motors and it is the decreasing density that is 
responsible for the monotonic decrease of J with w/;. 

In sharp contrast, at sufficiently low values of ujs/ujf, 
J varies non-monotonically with ujfi/ujf. In this regime 
of LUs/uf, at LUfi ^ 0, the high density of p causes steric 
hindrances which, in turn, leads to small J. When luji 
is "switched on" , p decreases with increasing lu fi and J 
increases up to a maximum because of the weakening 
of the hindrance effects. But, beyond a certain range 
of LUfi/iUf, the density of motors becomes so low that 
the movement of the motors is practically free of mu- 
tual hindrance; the decrease of J beyond its maximum is 
caused by the further reduction of density. Larger differ- 
ence between the predictions of our approximate analyt- 
ical calculations and computer simulation data at lower 
values of ujs/'-^f arises from the fact that the mean-flcld 
approximation neglects correlations which increases with 
increasing densitity of the motors. 

The above interpretation of trends of variations of J in 
Fig. [2] in terms of the corresponding variation of p is con- 
sistent with the results shown in Figs. [3] and [H But, why 
does p decrease monotonically with increasing ujfi/ujf? 
Increasing tofi, keeping all the other rate constants un- 
altered, leads to higher overall rate of transitions into 
strongly-bound states. Since, detachments of the motors 
from the microtubule track take place from the strongly 
bound state (see footnote [l^), the steady-state density 
is lower for higher values of w/j/w/. 

An approximate expression for J, which is obtained by 
retaining only the terms upto the first order in ujfi in a 
Taylor expansion of the right hand side of ([5]), is given 

by 
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J = Jo 



2Vd{1 + K] 



[(1 + K)ns + Kf + [(1 + K)n, - VDf - [iih - v^f - [n,, -Kf]+ o(w^j)(9) 



where Jq is the flux corresponding to ujfi — (i.e., in the 
absence of lane changing) . The approximate formula ([9]) 
still provides a reasonably good estimate of the flux even 
when ujfi is as large as ujf. 

In this communication we have extended the NOSC 
model for KIFl A traffic on MT [l^, by incorporating 
processes which correspond to shifting of the motors from 
one protofilament to another. These processes are analo- 
gous to lane changing of vehicles on multi-lane highways. 
On the basis of analytical treatment and computer simu- 
lations of the extended NOSC model, we have predicted 
the effects of such lane-changing on J, the steady-state 
flux of the KIFlA motors per lane. Over a wide region 
of parameter space, J decreases monotonically with in- 
creasing value of u!fi, a rate constant for lane-changing. 
However, in some regions of parameter space, J varies 
non-monotonically with increasing w/j. We have inter- 
preted the results by correlating the observed trends of 
variation of J with the corresponding variation of p, the 
average density of motors on a lane, and establishing the 



dependence of p on ujfi. 

Double-headed conventional kinesin rarely changes 
Double-headed dyneins may change lane in- 
|26| but, perhaps, not in-vivo [23|. The bound 
head of a double-headed motor imposes constraints on 
the stepping of the unbound head. Since such constraints 
do not exist for single-headed kinesins, KIFlA may find 
it easier to change lane. However, KIFlA may dimerize 
in-vivo [28j . Therefore, in-vitro experiments with fluo- 
rescently labelled KIFlA would be able to test our theo- 
retical predictions. In particular, variations of J with K 
and ujh (see figs El and S]) can be probed by varying con- 
centration of KIFlA and ATP molecules, respectively, in 
the solution. 
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